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We report the observation of prompt J/tp via double cc production from the e + e~ continuum. In 
this process one cc pair fragments into a J/ip meson while the remaining pair either produces a bound 
charmonium state or fragments into open charm. Both cases have been observed: the first by study- 
ing the mass spectrum of the system recoiling against the J /ip, and the second by reconstructing 
the J/tp together with a charmed meson. We find cross-sections of a(e + e~ — > J/ip r/ c (j)) x B(r\ c — > 
> 4charged) = (0.033 tg;^ ± 0.009) pb and a{e + e- -> J/ipD* + X) = (0.53 ± 0.14) pb, and 
infer a(e+ e - -> J/ip cc)/a(e + e~ -> J/tpX) = 0.59 + H ± 0.12; in each case the uncertainty quoted 
second is the systematic error. These results are obtained from a 46.2 fb _1 data sample collected 
near the T(4S) resonance, with the Belle detector at the KEKB asymmetric energy e + e~ collider. 



PACS numbers: 13.65.-K, 13.25.Gv, 14.40.Gx 
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Prompt charmonium production in e + e~ annihilation provides an opportunity to study both perturbative and non- 
perturbative effects in QCD. In a previous paper Q] we presented cross-section measurements for prompt J/ip and 
ip(2S) production, and studies of their kinematic properties, which were compared to predictions of the NRQCD B 
||, [|, || model. The BaBar Collaboration has also published results on prompt J/ip production B. NRQCD predicts 
that prompt J/ip production at \fs i=a 10.6 GeV is dominated by e + e~ — * J/ijjgg, with additional contributions from 
J/ip g, J/ipcc and other processes. The color-octet J/ip g signal predicted in Refs |Q, || is not observed |l| . The results 
in Refs (TJ [| do not constrain the contributions from J/ip gg and J/ip cc. 

We present the results of a search for J/ipcc production, i.e. double cc production, where the additional cc pair 
fragments into either charmonium or charmed hadrons. Exclusive associated production with a J/ip meson was 
proposed in Ref. as a method for searching for the remaining undiscovered charmonium states. Based on a non- 
perturbative approach the model of Ref. Q predicts a cross-section of a few pb for J / ip rj c production at y/s « 10.6 GeV. 
In contrast, recent theoretical estimates of the J/ipcc cross-section are as small as 0.07 pb |3|, || |9j. 

This analysis is based on 41.8 fb -1 of data at the T{AS) and 4.4fb _1 at an energy 60MeV below the resonance, 
collected with the Belle detector at the KEKB asymmetric energy storage rings jl0|]. Belle is a large solid-angle 
magnetic spectrometer jy]]. Charged particles are reconstructed in a 50-layer central drift chamber (CDC) and their 
impact parameters are determined using a three-layer silicon vertex detector (SVD). Kaon/pion separation is based 
on specific ionization (dE/dx) in the CDC, time of flight measurements and the response of aerogel Cerenkov counters 
(ACC). Electron identification is based on a combination oidE/dx measurements, the ACC response and information 
about the shape, energy deposit and position of the associated shower in the electromagnetic calorimeter (ECL). Muon 
identification is provided by 14 layers of 4.7 cm thick iron plates interleaved with resistive plate counters. Photons are 
reconstructed in the ECL as showers with an energy larger than 20 MeV that are not associated with charged tracks. 

We use hadronic events separated from QED, rr, two-photon and beam-gas interaction backgrounds by the selection 
criteria described in Ref. [Q . Charged pion and kaon candidates are well reconstructed tracks that have been positively 
identified. Ks candidates are reconstructed by combining ir + ir~ pairs with an invariant mass within 10MeV/c 2 of 
the nominal Ks mass. We require the distance between the pion tracks at the Ks vertex to be less than 1 cm, 
the transverse flight distance from the interaction point to be greater than 0.5 cm and the angle between the Ks 
momentum direction and decay path to be smaller than 0.1 rad. Photons of energy greater than 50 MeV are combined 
to form 7T° — > 77 candidates if their invariant mass lies within 10MeV/c 2 of the nominal ir° mass. Such 77 pairs are 
fitted with a ir° mass constraint to improve the momentum resolution. 

The J/ip — » £ + £~ reconstruction procedure is identical to that presented in Ref. Q. Two positively identified 
lepton candidates are required to form a common vertex that is less than 500 /an from the interaction point in 
the plane perpendicular to the beam axis. For J/ijj — > e + e~ , the invariant mass calculation includes the four- 
momentum of photons detected within 50 mrad of the directions, as a partial correction for final state radiation 
and bremsstrahlung energy loss. QED processes are suppressed by requiring the total charged multiplicity (N<±) to 
be greater than four. J/ip mesons from BB events are removed by requiring the momentum of the meson in the e + e _ 
center of mass (CM) system, p}/^, to be greater than 2.0GeV/c. The process e + e~ — > ip(2S)~f is partially removed 
by rejecting events with a photon of CM energy E* > 3.5 GeV. 

The J/ip — > £ + £~~ signal region is defined by the mass window \Mg+ t - — Mj/^\ < 30MeV/c 2 (« 2.5cr). The 
sideband region is defined by 100 < \M e + e - — Mj/^\ < 400MeV/c 2 and is used to estimate the contribution from the 
dilcpton combinatorial background under the J/ip peak. The sideband is scaled, parameterizing the background by a 
quadratic function. 

Distributions of the mass of the system recoiling against the J/ip candidate — the "recoil mass" — after this selection 
are shown in Fig. [j]a, for both signal and sideband regions. The recoil mass is defined as 



M rocoil = ^/(yi j - P *? r (i) 

A clear threshold near 2 m c can be seen. To avoid model dependence, we do not perform an acceptance correction. 
As examples, the Monte Carlo efficiencies for both e + e~ — > J/ipqq (q — u,d,s) and J/ipcc events are shown in the 
inset: the variation with M reco ;i is smooth. For recoil masses smaller than 2.8GeV/c 2 , no significant J/ip signal is 
observed: a fit to the dilepton mass spectrum for events in this region gives 16.4 ± 7.8 events. 

The region 2 m c < M rcco ji < 2 mu is studied in more detail in order to search for production of J/ip together with an 
additional charmonium state. We apply a mass constrained fit to the J/ip candidates before determining P*j^ ■ Monte 
Carlo simulation predicts that this improves the M rcco ji resolution by a factor close to two. The resulting recoil mass 
spectrum in the data is presented in Fig. ^p: a clear peak is observed around 3 GeV/c 2 . Since e + e _ — > 7* — * J/ip J/ip 
is forbidden by charge conjugation symmetry, we interpret this peak as e + e _ — ► J/ipr] c . Additional peaks at recoil 
masses consistent with the Xco and rj c (2S) mass |L3 are also seen. 
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FIG. 1: a) The recoil mass distribution for the J/ip signal region (points) and scaled sidebands (shaded histogram). The inset 
shows the MC reconstruction efficiency for e + e~ — » J/ip qq (q = u, d, s; open circles) and e + e~ — > J/ip cc events (closed circles), 
b) The recoil mass distribution after refitting the J/ip candidate with a mass constraint for the J/ip signal region (points) and 
scaled sidebands (shaded histogram). The curve represents the fit described in the text. 



To reproduce the rj c shape in the recoil mass spectrum we generate e + e~ — > 7* (7) — > J/iprj c (^) Monte Carlo 
events. We assume p-wave J/ip rj c production, as required by parity conservation, and phase space suppression of 
7* — > J/ipijc as the virtual 7* energy varies due to initial state radiation (ISR). This effect produces a high-mass 
tail in M rcco ji. The corresponding r/ c (2S) lineshape is narrower, due to the larger mass of the state and its presumed 
smaller intrinsic width. The Xco lineshape assumes s-wave production. The effect of varying the energy dependence 
of the cross-section in the ISR calculation is included in the systematic error. 

A fit to the recoil mass spectrum finds an r\ c yield N 71c = 671^2 & t a ma ss M = (2.962 ± 0.013) GeV/c 2 , with 
Xcd,r} c {2S) yields iV Xc0 = 39±^, N Vc(2S ) = 42±J| at masses M Xe0 = (3.403 ± 0.014) GeV/c 2 and M Vc(2S ) = (3.622 ± 
0.012) GeV/c 2 respectively. We use a second order polynomial to describe the background. Only the region below 
the open charm threshold (M reC oii < 3.73 GeV/c 2 ) is included in the fit. 

We assess the significance of each signal i using er^ = y/ —2 \n(C l / C max ), where C max is the maximum likelihood 
returned by the fit, and C l Q is the likelihood with the yield of the state i (i = rj c ,x c o,rj c (2S)) set to zero. We find 
= 6.7, <t Xc0 = 3.3 and (J T]c (2S) — 3.4. As the significance of the Xco an d rj c (2S) peaks is low, we perform an 
additional fit using only the r\ c shape and the polynomial background: this finds N Vc = 56^^ with c, )c = 5.9. We 
use this result, together with the results of fits after varying the charmonium intrinsic widths |l3| and the choice of 
background shape, to estimate the systematic error on the r\ c yield due to the fitting procedure. 

To determine the e + e~ — ► J/iprjc^) cross-section we correct the signal yield for the reconstruction efficiency 
obtained from the Monte Carlo. Because of the requirement iV c h > 4, the recoil system must contain at least three 
charged tracks: this removes r\ c decays into or 2 charged tracks plus neutrals. As r\ c branching fractions are poorly 
known, we express our result in terms of the product a(e + e~ — > J/ipT] c (7)) x B(r\ c — > > 4 charged), which we find 
to be (0.033 j^Q QQg ± 0.009) pb. Here and elsewhere, the uncertainty quoted second is the systematic error. The 
various sources of systematic error are listed in Table |, with dominant contributions from the uncertainty in the ISR 
correction, the choice of fitting procedure, and the production and helicity angle distributions for the J/ip meson. 

To study the J/ipcc mechanism in the region M reC oii > 2mfl, we search for fully reconstructed D* + and D° 
decays [jL4j in events with a J/ip meson. For the study of J / ip D* + associated production we reconstruct D* + — > D° n + 
using five D° decay modes: K~ir + , K~K + , K^tt^tt + tt + , Kstt + it~ and K~Tr + Tr°. We select D° candidates in a 
±10 MeV/c 2 mass window for the charged modes and a ±20 MeV/c 2 window for K~ir + ir° (approximately 2 a in each 
case). To improve the M D a^+ resolution D° candidates are refitted to the nominal D° mass. 

Although B — * J/ipX decays are rejected by the selection, semileptonic B decays contribute to the background 
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TABLE I: Sources of systematic error for a(e + e — > J/ipr) c ). 



Source 


Systematic error (%) 


ISR correction 


±19 


Fitting procedure 


±16 


J/ip polarization 


±11 


Track reconstruction 


±5 


Lepton identification 


±4 


Total 


±28 



under the J/ip peak and lead to a large D* + signal. To remove the remaining BB background we require that either 
the D* + or one of the leptons from the J/ip candidate have a momentum above the kinematic limit for B decays: 
p* D ,+ > 2.6GeV/c or p* (± > 2.6GeV/c. Choosing the D°ir + combination with the best D° mass yields at most one 
J/ipD* + candidate per event. 

The scatter plot of the dilepton mass versus the D°w + mass, and the D°ir + mass projection, are shown in Figs. |]a,&. 
We perform a fit to the D°ir + mass distribution in the J/ip signal window, with a Gaussian for the D* + signal and 
a threshold function A^/AI D a 7I + — A/ t hrcs for the background. The J/ip sideband is fit simultaneously, and used to 
estimate the combinatorial (£ + £~) D* + contribution to the D* + yield in the signal region. We find N D *+ = 10.5^' q, 
with a combinatorial contribution of 0.4 ± 0.3: the signal yield is lO.ltg q, with significance CTj/,/, £>*+ x — 5.3. As a 
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cross-check we fit the dilepton mass distribution in the region 2.008 < Mno n + < 2.012 GeV/c 2 , finding Nj/^ = 9.6_ 2 9 
The J/ip signal shape is fixed from the Monte Carlo simulation and we include a linear background function. 
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FIG. 2: Results of a search for associated production of J/ip and charm mesons: a) the scatter plot M(l + l ) vs M(D°7r + ); b) pro- 
jection onto the Al(D°n + ) axis; c) the scatter plot M(l + T) vs M(K~n + (K + K~)); d) projection onto the M(K ~n + (K + K~)) 
axis. Points with error bars show the J/ij) signal region and the shaded histograms show the scaled sidebands. The curves 
represent the fit described in the text. 



For the study of J /ip D° associated production we use only the cleanest D° decay modes D° — > K~ir + and K~K + . 
As in the J/ip D* + study, we remove BB events by requiring p* D „ > 2.6 GeV/c or p* t± > 2.6 GeV/c. 

A plot of dilepton versus K~ir + (K~K + ) masses, and the projection onto the K~ir + (K~ K + ) mass axis, are shown 
in Figs. ||c,rf. A simultaneous fit to the K~n + (K~ K + ) mass distribution in the J/ip signal window and the sideband 
finds Npo = 15.9~^; 7 in the signal region, with a combinatorial (£ + £~) D° contribution of 1.0 ± 0.8. The signal yield 
is 14.91"!' 8' w ith significance (Jj/^,d°x = 3.7. We use a Gaussian signal shape and a linear background function 
in the fit. As a cross-check, we also fit the dilepton mass spectrum for D° signal and sideband regions, obtaining 
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Nj/ip = 17.7l 4 ; 6 and Nj/^ = 4.3 ± 0.8 respectively, where the sideband number is the result of the fit scaled to the 
expected contribution under the D° peak. For all fits the signal shapes are fixed from the Monte Carlo simulation. 

To study our reconstruction efficiency for J/ip mesons produced together with a charmed meson, we generate 
e + e" — > J/ipcc Monte Carlo events using a simple model adapted from the QQ event generator Since the 

efficiency for both particle reconstruction and selection criteria strongly depend on the kinematics, we correct the 
general kinematic characteristics of the Monte Carlo events to match those of the data using a large sample of 
continuum J/ip events. In particular, the distributions of the recoil mass, the J/ip production and helicity angles, and 
the angle between the thrust axis of the cc system and its boost are adjusted to match the data. These quantities 
almost fully describe the kinematics of e + e" — * J/ipcc. The fragmentation function of cc into charmed mesons in 
each bin of Q 2 (cc) is the only characteristic that has an effect on the efficiency that cannot be determined from the 
data. We vary the fragmentation function over a wide range, taking the difference in the efficiency into account as a 
systematic error. 

The efficiency is first calculated for Pj/^, > 2.0 GeV/c and then extrapolated to the full momentum interval, taking 
into account the cross-sections for inclusive continuum production obtained in Ref. Q for both Pjm > 2.0 and p*j ,^ < 
2.0GeV/c data. The overall efficiencies for J/ipD* + and J/ipD° are calculated to be 6j/^, D *+ = (4.1 ± 1.0) ■ 10" 4 
and ejupo = (3.7 ± 0.8) ■ 10" 4 respectively. Using these values, we find cross-sections er(e + e~ — > J/tp D* + X) = 
(Q-53 ±8:ii ± 0.14) pb and cr(e+e" ~> J/ipD°X) = (0.87 to.ll ± 0.20) pb. Contributions to the systematic error are 
summarized in Table [H]. 

According to the Lund model, cc fragmentation produces charmed mesons at the rate of 0.53 per event for D* + , and 
1.18 per event for D°, where both numbers include feed-down from higher states (in particular, D* + — > D°tt + ) [[hs). 
Assuming that these rates apply to cc fragmentation in e + e~ — * J/ipcc, we calculate cr(e + e~ — > J/ipcc) and find 
(1.01 toil ± 0.26) pb and (0.74 toil ± 0.19) pb based on our D* + and D° measurements, respectively. No systematic 
error is included for our use of the Lund fragmentation rates. These results are slightly correlated, as two events 
are common to both samples. Taking this into account, we average the results and obtain (0.87 Iq 19 ^ 0.17) pb. In 
Ref. |l[ we found the inclusive prompt J/ip cross-section to be <r(e + e~ — > J/ipX) = (1.47 ± 0.10 ± 0.11) pb, based 
on a 32.71b" 1 dataset. We therefore infer that a large fraction of prompt J/ip events, cr(e + e~ — > J/ipcc)/a(e + e~ — > 
J/ipX) = 0.59 

-o'.i3 ^ 0-12, is due to the e + e — > J/ipcc process. Contributions to the systematic error on this ratio 
are shown in Table [n| by the numbers in parentheses. 

This J/ipcc cross-section is an order of magnitude larger than predicted in Rcfs. ||, || and contradicts the 
NRQCD expectation that the J/ipcc fraction is small ||, ||]. We note, however, that our result is dependent on 
the fragmentation model assumed for the cc system. In the future, more comprehensive measurements including 
associated J/ipD + , J/ipDf and J/ipA+ production could significantly reduce this model dependence. 

TABLE II: Systematic error contributions for a(e + e~ — » J/ipD X). Numbers in parentheses show contributions to the error 
on the ratio a(e + e~ — ^ J/ip cc)/a(e + e~ — > J/ipX). 





Systematic error (%) 


Source 


J/ip D° 


J/ip D*+ 


MC kinematics correction 


±11 (±8) 


±10 (±8) 


cc fragmentation function 


±8 (±8) 


±15 (±15) 


Fitting procedure 


±10 (±10) 


±5 (±5) 


Efficiency of pj/^, cut 


±11 (0) 


±11 (0) 


Track reconstruction 


±8 (±4) 


±12 (±8) 


Lepton and K identification 


±6 (±3) 


±6 (±3) 


Total 


23 (16) 


26 (20) 



In summary, we have observed both a charmonium state and charmed mesons accompanying prompt J/ip production 
in e+e" annihilation. We measure cr(e + e" -> J/tpr) c (j)) x B{r/ C -> > 4 charged) = (0.033 to.ool ± 0.009) pb and 
CT(e+e" J/ipD* + X) = (0.53 ± 0.14) pb, and estimate o-(e+e" -> J/ip cc)/a{e+ e - -► J/ipX) = 0.59±g;i| ± 
0.12. Our results favor e + e" — * J/ipcc as the leading mechanism for prompt J/ip production at i/s w 10.6 GeV. 
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